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Abstract. The indirect RKKY interaction in iron pnictide and chalcogenide metals 
is calculated for a simplified four bands Fermi surface (FS) model. We investigate 
the specific multi-band features and show that distinct length scales of the RKKY 
oscillations appear. For the regular lattice of the local moments, the generalized RKKY 
interaction is defined in momentum space. We consider its momentum dependence in 
paramagnetic and spin density wave (SDW) phases, discuss its implications for the 
possible type of magnetic order and compare it to the results obtained from more 
realistic tight-binding type Fermi surface model. Our finding can give important clues 
on the magnetic ordering of the 4f- iron based superconductors. 
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1. Introduction 

The discovery of iron-based superconductors [1], has lead to a renewed interest in 
multiband superconductivity [21 [3] . Their parent compounds are paramagnetic metals 
at high temperature and mostly show an antiferromagnetic (AF) order of spin density 
wave (SDW) type slightly below a tetragonal to orthorhombic structural transition j4]. 
Magnetic order of the itinerant Fe moments is driven by the nesting properties of 3d type 
Fermi surfaces which consist of hole and electron pockets around the center and at the 
boundaries of the Brillouin zone (BZ). The magnetism of the parent compounds becomes 
even more involved in rare earth (R) based Fe-pnictides where layers with localized 4f 
moments exist that are separated from the Fe-As layers with itinerant 3d moments. It 
is found that the R- moments order at temperatures much below the SDW transition. 
Their magnetic structure may be generally different from that of the Fe layers. The 
latter predominantely show the collinear in-plane stripe structure of Fe moments while 
the R layers may order in different noncollinear or simple ferromagnetic structure. For 
example the in-plane FM order appears in the 122 family RFe 2 As 2 for R = Eu[5j, and 
in the 1111 family RFeAsO the R = Ce compound shows non-collinear in-plane stripe 
structure [6l [7j. 

Because of the close connection of magnetic and electronic properties in these kind 
of materials, investigation of magnetic order and effective coupling mechanism of local 
moments is a question of great interest. The latter may belong to a periodic sublattice of 
4f-elements [HI El El [5] ; doped disordered impurity magnetic moments (diluted magnetic 
moments) [10l [HI [12], or possibly partially localized 3d electrons on the Fe sublattice 
originating from an orbitally selective Mott-Hubbard localization. 

A well studied example of the 4f- based 122 Fe pnictides is EuFe 2 As 2 , with its 
highest SDW transition temperature reported at T SDW = 190K [51 H51 fT4l H51 [T6l H7| . 
The AF ordering of localized 4f- moments (Eu 2+ with S = 7/2) happens at lower 
temperature, T N « 20K, for a different wave vector Q = (0, 0, 1). This means the Eu 2+ - 
moments are FM ordered in a6-planes which is in contrast to the columnar AF order of 
Fe itinerant moments [5j. In addition, the relaxational behavior of Eu 2+ spins in ESR 
shows a distinct magnetic anisotropy below Tsdw |1SJ . The magnetization anisotropy 
of Eu spins has a temperature dependency which is changing across the SDW transition 
temperature [12]. These observations indicate that the SDW transition in the Fe-3d 
itinerant subsystem may have important effects on the effective coupling of local 4f 
moments. Evidence for the coupling of the 4/ 7 -electron and conduction electrons is also 
obtained by ESR in the Gd-based 1111 compounds La x Gdi_ x FeAsO (Gd 3+ with spin 
S = 7/2) [12]. Below T ~ 6K, neutron scattering measurements also found the Fe-Nd 
interaction in NdFeAsO single crystals which force the Nd- moments to be ordered AF 
with the same a6-plane configuration as Fe moments [20]. 

Since the 4f moments are localized, their direct exchange can be neglected. How- 
ever one has to consider the indirect Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion [211 E21 [23] of 4f moments via spin polarization of 3d conduction electrons. The 
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Figure 1. (Color online) a) The Fermi surface with pure electron- like pockets with 
ellipticity e = 0.4, centered around the Q/3 [Qi = (±7r, 0) and Q2 = (0, ±7r)] points in 
the unfolded BZ (/3-bands). b) Fermi surface with pure hole- like pockets centered at 
T- point (a-bands). c) Fermi surface with two circular hole pockets and two electron- 
like pockets. Note that c) results from a combination of electron and hole pockets in 
a) and b). 



strength of this oscillatory effective exchange is controlled by the distance between two 
localized 4f moments and the Fermi surface (FS) properties of 3d conduction electrons. 
Such an interaction generally plays an important role in revealing the nature of the 
magnetism in metals with partially filled d- and /-electron shells. It is well-known that 
the RKKY interaction for an anisotropic Fermi surface (FS) with nesting consists of 
several terms originating from flat regions in the conduction bands (van Hove regions) 
and those which describe the interference between contributions from their vicinities 
[24] . For the nearly nested conduction bands the latter term is present down to the 
interatomic distances and favors the commensurate antiferromagnetic ordering of the 
localized moments. The degree of nesting varies continuously for different compounds 
[25] , where we have a perfect nesting in some pnictides [26] and weakly nested compound 
of some chalcogenides [271 [28] . 

This motivates us to investigate how the RKKY mechanism in iron based 
superconductors depends on the absence or presence of electron- or hole- Fermi surface 
sheets and their nesting properties. In most of the iron based superconducting 
systems, the Fermi surface consists of the hole and electron pockets simultaneously 
|2Hl EQl EH E21 [331 EH [26]. It is also possible to have a situation with only the 
hole- like or the electron- like pockets separately. For example, according to band 
structure calculations [35j and angle resolved photoemission spectroscopy (ARPES) [36] 
the overdoped KFe 2 As 2 compound has no electron - like FS sheets [37]. On the other 
hand ARPES studies on the alkali-intercalated Fe chalcogenide system KFe 2 Se 2 have 
recently shown that there are no hole- like FS sheets for some doping range [38]. We will 
use a simplified parabolic band model with isotropic hole but anisotropic electron mass 
and study how the real space oscillations of RKKY local moment pair interaction evolve 
with such FS change. We study the general case involving four bands that consist of two 




Figure 2. (Color online) The fully analytical calculations of the individual components 
of magnetic spin susceptibility of conduction electrons x 77 ( r ) (obtained using Eqs.([l2j 
fl3|) along the x-direction for a) pure electron- like pockets, b) pure hole- like pockets. 
They are completely in agreement with the numerical results, c) The numerical 
calculation of the individual components of magnetic spin susceptibility of conduction 
electrons x 77 ( r ) along the x-direction for only the inter hole- and electron- like 

/ / 2 

contributions. Note that x 77 ( r ) = X 7 7 ( r )? furthermore xo = m ^ • 



electron- like and two hole- like sheets. We compare these results with the special case 
where only one type of sheet is present and the closed analytical solution is obtained. 

For a lattice of local moments with effective RKKY interactions the latter will 
also be analyzed in momentum space. This is relevant for the R-based pnictides 
with magnetic order appearing both in the itinerant 3d and localized 4f subsystems. 
These results will be compared to a more detailed calculation using tight-binding type 
conduction electron bands. As indicated by the above examples the SDW transition and 
associated gap opening may influence the q- dependence and anisotropy of the effective 
RKKY exchange interaction which determines the type of order in the local moment 
subsystem. This feedback effect on local moment magnetism will also be investigated 
for the tight binding model. 

2. Model definition 

The Hamiltonian describing localized magnetic moments in the multi-band conduction 
electron sea is given by 

% = He + Hint, (1) 

where % c is the conduction electron Hamiltonian according to 

= £ k C \ka C 7ka' (2) 

k,7,cr 




Figure 3. The normalized total RKKY interaction, J n (r)/Jo, for a) pure electron- like 
pockets, b) pure hole- like pockets and c) for the case with 4 electron and hole bands, 
corresponding to the case of Fig. ([He), along the x-direction, furthermore Jo = JexXo- 
(a) and (c) exhibit rapid oscillations due to the inter-pocket contributions (c. f. 
Fig. EM 




Figure 4. (Color online) The spatial plot of the normalized total RKKY interaction, 
j n (r) = i j (r)|+o 015 •> f° r a ) P ure electron- like pockets, b) pure hole- like pockets 
and c) for the case with 4 electron and hole bands, corresponding to Fig.(pQ). 
Bar plots: (d), (e) and (f) present of normalized total RKKY interaction, (j n (r)), 
for two local moment sitting on the origin and the lattice site (x, y) = (na, ma); n, m = 
0,=bl.. with distance d — \J x 2 -\-y 2 corresponding to the plots (a), (b) and (c), 
respectively. 
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Here, c^ kcr denotes the creation operators and is the dispersion of conduction electrons 
with momentum k = (k X) k y ) and spin a in band 7. 

In the present work we consider primarily the minimal four-bands model consisting 
two hole- like pockets centered at T- point (a-bands), and two elliptical electron FS 
pockets centered around [Qi = (±7r, 0) and Q 2 = (0, ±7r)] points in the unfolded 
Brillouin zone (BZ) (/3-bands). We consider three different cases: a) only electron- like 
FS sheets, b) only hole- like FS sheets, and finally c) both hole and electron FS sheets 
(see FigH]). Thus one can rewrite the total conduction electron Hamiltonian Eq.([2]) as 

= Y £ k a ak<T a aka + Y ^k^ko-^ka, (3) 
a,k,cr /5,k,cr 

where a* k<7 (fel^) creates electrons with spin a, momentum k at hole- (electron-) like 
band a (/?). The dispersion of the electrons pockets can be modeled as ^k-Q^ = 
7^^ + — — a ^ for the hole- like pockets we have = — [^-^ + ^-^-] — /x a , where 
denotes the band mass [n = x, y) in k n direction. Furthermore we define \i$ = \± + £ , 
and fi a = /a — £0 where s and \± are the energy offset and the chemical potential, 
respectively. We consider for the hole- like pockets = m a with m a2 = 3m ai = 3m e , 
and for the electron- like pockets m^ 1 rriy 2 — (1 + e)m e , m^ 1 m^, 2 = (1 — e)m e with 
ellipticity e and electron mass m e . The corresponding general FS is shown in Fig. ([He). 

The interaction between the spin of conduction electrons, s(R), and the moment 
of localized /-electrons, S^, at site is given by 

U int = J ex £ J(R - R,) s(R) • S u (4) 

where the exchange integral /(R — R^) is established by the overlap of 3d conduction 
and localized 4f electron wave functions. The real overlap of the conduction electron and 
localized 4f electron wave- functions and therefore /(R — R^) in general is non-zero in a 
finite volume around the 4f site R^. This leads to a damping in the effective intersite 
coupling energies [391 [40l [41] . For simplicity we use here the common approximation of 
an on-site exchange interaction which can be approximated as /(R— R^) = J ex S(H— R^). 
Here J ex is the on-site exchange coupling constant that can be obtained from a more 
microscopic Anderson-type model [42] by a Schrieffer- Wolff transformation. 



3. Spatial and momentum dependence of effective RKKY interaction 

The contact exchange interaction with local moments will polarize the conduction 
states which then leads to an effective RKKY exchange between moments at different 
sites. It may be obtained by using the standard second-order perturbation theory with 
respect to Hint- Here we consider first the general form of the RKKY interaction of 
a pair of local moments and its distance dependence for the various FS topologies. In 
addition we investigate the momentum dependence of the effective RKKY exchange 
which determines the magnetic structure in a periodic lattice of local moments. 
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3.1. Localized magnetic moment pairs and the distance dependence of RKKY 
interaction 

In the paramagnetic or normal state regime the effective exchange Hamiltonian of the 
RKKY interaction which describes the interaction between two local impurity spins at 
the positions i and j is derived as (2H 143] 

Wrkky = —Jn(?)Si ' Sj, (5) 
where r = — Rj = xx + yy ) and the effective exchange couplings are then given by 

Mr) = £ J^' (r) = J 2 ex Re ( £ e^-^ X ^\v) J , (6) 
7 y y 7 y y 

here x 11 ' { Y ) denotes the intra- (7 = 7') and inter- (7 7^ 7') magnetic spin susceptibility 
of conduction electrons (Lindhard response function) which is defined by 

X 77 '(r) = -k B T J2 C 7 (r, iu n )Gy (r, iu n ). (7) 

n 

Here uj n = nT(2n + 1) is the fermionic Matsubara frequency and 

G 7 (k, iu n ) = - — - — ? , (8) 

is the conduction electrons Greens function in momentum representation. After some 
algebra [44J the real space Greens function can be obtained as 



rnirriy 



G^Y.iUn) = 2 if (^2Z 7 (a; n )p 7 ), (9) 

where K (. . .) is the modified Bessel (Macdonald) function. Here we define Z 1 = ji 1 -\-iuj n 
and p 1 = (m],x 2 + m^y 2 )/h 2 . In the low temperature regime 



rioo 

**rE(...)->— /. M-), 

m LTX I J —too 



therefore the magnetic spin susceptibility, Eq.([7]), can be calculated using the Green's 
function as [15] 

1 r io ° 

X 11 {r) = — duG^r, jw)Gy(r, iw). (10) 

J —ioo 

The closed analytic expression for the above integral may be obtained in the case 
of having the same kind of bands, where /x 7 = fly (FigUla or FigHJb). Explicitly, using 
the properties of the modified Bessel functions in the limit of Z 1 —> /x 7 + zO ± , we have 

K n (j2Z^) = ^e^Ht\^MhP,\ (11) 



where H^\x) = 3 n (x) ± iY n {x) are Hankel functions 

In the special case p 1 = py, by using Eqs. (l9l[T0l) . the magnetic susceptibility can 
be written as 



X 77 '(r) = $ 



77' 



jotx^Yopg + j^x^pg 



(12) 
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(a) (b) 




(c) (d) 

Figure 5. (Color online) The momentum dependency of the Fourier transformation 
of the total RKKY interaction, J n (q) for a) pure electron- like pockets, b) pure hole- 
like pockets c) only the inter hole- and electron- like contributions and d) the to total 
interaction corresponding to sum of all FS sheet contributions. Results are obtained 
for the simple parabolic dispersion corresponding to Fig.(pQ). Absolute maximum of 
J nisi) m this simplified model occurs at q = and side maxima at q = Q a . 



here $ 7 y 



mZmyml my ^4 and X 1 J\2p 7 p 7 \. 



In the general case p 1 ^ py we still can find the closed solution 

X 77' (r) 



^77 



y 7 2 - y 7 2 , 



x 



- Y 7 /(j!(Y 7 /)Y (Y 7 ) + Yi(Xy)J (X 1 



(13) 



For the simple case of spherical FS sheets, i.e. m], = rri^ we have X 1 = kjr, where 
kj is the Fermi momentum, and r = ^Jx 1 + y 2 the distance between local moments. 
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Figure 6. (Color online) (a) A typical tight binding band structure for iron based 122- 
type superconductors with two circular hole pockets and two electron- like pockets. The 
BZ symmetry points are T(0,0), X(0, f ), F(f , 0) and M(£, f ). The band structure 
parametrization can be found in Ref.[47j. (b) The corresponding density of states 
(DOS), and (c) its Fermi surface, (d) The reconstructed FS in the SDW phase with 
ordering parameter W = 0.15ey. 



3.2. Periodic lattice of local moments and the momentum dependence of RKKY 
interaction 

In the case of local moments forming a regular lattice it is necessary to go to the 
momentum representation of the RKKY interaction which is defined by 

«// = - E U*)*i • Sj = - £ </n(q)S_ q • S q , (14) 

ij q 

where the real space RKKY interaction is given by 

Ur) = ^E e ' (Q7 " Qy) ' r X 77 '(r), (15) 

and 

Uq) = </ e lEx 77 '(q- (Q7 - Qv)), (is) 

is its Fourier transform and S q are the Fourier components of S^. They are defined by 
U<D = ^= E Jn(r)e^ r , <? q = -^=£ Ste** (17) 
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respectively, with TV denoting the lattice size. The RKKY interaction is proportional 
to the static conduction electron magnetic susceptibility given by 

iV k ^k-q - ^k 

where /(...) denotes the Fermi function. Here we implied that the electron pocket 
has been shifted to the T- center of the BZ. We may also conveniently define a total 
susceptibility 

x t (q) = Ex 77 '(q-(Q7-Q7'))- (19) 

The meaning of x*(q) is twofold: At the wave vector q where Xt(q) acquires its maximum 
an SDW instability may occur due to 3d conduction electron interactions, opening a 
SDW gap at FS patches connected by this momentum. Furthermore at the same wave 
vector the ordering of 4f local moments due to J n (q) = Je X Xt(<l) i n the rare earth - 
based Fe pnictides should be expected at lower temperature. Note, however, that due 
to the feedback effect the opening of the SDW gap may influence the wave vector where 
J n (q) has its maximum. 



3.3. Momentum dependence of RKKY interaction in spin density wave phase 



Recently it was shown [43] that the change of the FS nesting (resulting of the opening 
of the SDW gap) influences both strength and oscillatory behavior of RKKY in the 
antiferromagnetic state. In the mean-field approximation the SDW ordering can be 
described by [48] 



Hsdw = E W ° hLw Ql . + h.c 



(20) 



ka 



where W is the ordering amplitude and the spin index a — ± corresponds to spin f and 
^. The SDW ordering induces an anisotropy in the RKKY interaction and maps the 
RKKY into an effective anisotropic XXZ-type Heisenberg exchange model, 



H 



ff 



E 
E 



q l 



J x (r)(s?S? + S?S?)+J a (r)S?S 
J,(q)(S*S*_ q + S»S» q ) + J,(q)S q S 



(21) 



where the momentum dependence of the Fourier transform of the total RKKY 
interaction in SDW phase is obtained as 



j (a) - /2 V ^ z ^- q )-/(£k) 



77 x ,k 



^k-q 



(22) 



with the quasiparticle energies defined by El = e£, E£ = e£ 2 , and 



E, 



1,2 



(4 1 + 4 1 )±V(4 1 -4 1 ) 2 + 4^ 
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(c) (d) 

Figure 7. (Color online) The momentum dependence of the Fourier transform of 
the total RKKY interaction, J n (q), for a) pure electron- like pockets, b) pure hole- 
like pockets, c) only the inter hole- and electron- like contributions and d) the total 
interaction corresponding to sum of all contributions. Results correspond to a typical 
tight-binding type band structure example (see Fig.flBJ)) for 122 compounds. 



The coherence factors ^y kq are given by 



77 x ,k,q 



1 k-q,+ 1 k,- 



and 



77 x ,k,q k— q,+ k,+ ' k— q,- 

_ >w 



\b z - T 7 T 7 

with T^'J = (u k ± av k ) 2 and Tj^ 
unitary transformation are given by 



Ti 



T 7 
i k,-5 



(23) 



(24) 



2 2 



l ± 



1. Furthermore the coefficients of the 



(25) 



4. Discussion of numerical results 



Now we begin our numerical discussion of the real space and momentum dependence 
of the RKKY interaction. For the calculations we use /x 7 = h 2 kj 2 /2m e with typical 



Evolution of the multiband RKKY interaction 



12 



values of Fermi momentum kja = O.lbir for hole (a) - like pockets, and k^a = 0.2n with 
ellipticity e = 0.4 for electron (/?)- like pockets (FigHJ). Here a is the lattice constant. 

First we present the results for the individual intra- and inter-band contributions 
to the total spin susceptibility for the various FS models in Fig.(pQ), they are shown 
in Fig. (EI) in the same sequence. The distance dependence of individual components 
of x 77 (r), along the x-direction is shown in this figure. We notice that the oscillation 
pattern is similar for intra-band contributions (a — a, (3 — /3), however their maxima and 
minima are slightly shifted due to the different FS dimensions visible in Fig.(pp). This 
is particularly the case in the inter-band contributions of Fig.([2lc). We note that the 
results in Figs.([2la,b) were obtained from the closed analytical expression in Eq. f fT3l) 
which are identical to the numerical results. For the general FS case with both pockets 
present x 77 ( r ) can only be calculated numerically. 

The total RKKY exchange is proportional to the sum of all (77') contributions 
multiplied with a phase factor determined by the pocket distance, i.e., the nesting 
vectors Q 7 — Qy. It is shown in Fig. (ED for the three Fermi surface cases of Fig.(pp). 
Fig.(|3la) represents the RKKY interaction for pure electron- like structure. As a result 
of the extra inter-band (13-/3) phase factors in Eq. f |T5l ). it shows a new rapid superposed 
oscillation on top of the fundamental oscillation determined by the FS sheet diameters 
which is created by summation of individual parts of Fig.([2la). 

For the second case with only hole- like pockets, the total RKKY coupling shows 
only the oscillations defined by total summation of individual terms. In contrast to 
the case of electron pockets the rapid oscillations are absent because both (a 1? a 2 ) hole 
pockets are T- centered without a shift between them. 

Finally we present the calculated RKKY interaction for general four band case 
along the x-direction in the Fig.([3lc). In similar way than in the pure electron case, 
it shows the combination of slow overall and additional rapid oscillation. The latter 
originate from both electron-electron ((3 — (3) and electron-hole (a — (3) contribution but 
not from the hole-hole (a — a) part. The size of the electronic pockets and the corre- 
sponding Fermi momentum can be changed by tuning the chemical potential. Since the 
oscillation is defined by the Xy(— kjr) value, an increase of kj causes the reduction of 
the wave length of both oscillation types and vice versa. 

For a better understanding of the RKKY oscillatory behaviour, in Fig.(|Ha-c) we 
show the full spatial dependency of the normalized total RKKY interaction, j n (r) = 
\j J (r)\+r] ( we use P arame ter r] = 0.015 for good extremal contrast). Fig.(|Ha) is 
for pure electron- like pockets, Fig.((4lb) is for pure hole- like pockets and Fig. ((He) is 
for the general case with 4 electron and hole bands. For a FS with only hole pockets 
(Fig.(|4lb) only long range (radial) oscillations which are almost isotropic appear. When 
electron pockets are present (Fig.(@ia,c) there are superposed short range oscillations 
(also azimuthal) due to the interband processes. However a shell-like overall structure 
of FM/AF regions for moderate distances is preserved, with a notable phase shift by a 
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half period along the (1,1) direction. 

One of the important issues of the diluted system with random local moments (e.g. 
magnetic impurities like Mn) sitting on lattice sites is the conditon for their magnetic 
ordering. For this purpose it is useful to know the FM/AF oscillatory RKKY behaviour. 
Therefore we present the normalized total RKKY interaction, j n (r), for two local mo- 
ments sitting at the origin and the lattice site (x,y) = {na^ma)\n^m = 0,±1... with 
a distance d = yjx 1 + y 2 in the plots of Fig.(|Hd-f), respectively. These plots clearly 
show the change of the interaction form AF to FM by varying the distance of the local 
moment impurities which is much more rapid when electron like pockets are present. 
The average impurity distance can be controlled by the concentration of the magnetic 
impurities in the sample [TT]. 

An experimental determination of the real space variation of the RKKY interaction 
is sofar not easily possible. However, its Fourier transform is a more accessible quan- 
tity, because the magnetic ordering both in the 3d system and in the effectively RKKY 
coupled 4f system should occur at wave vectors that are maxima of the total static sus- 
ceptibility x*(q) = Jn{o)/ Jl X ) therefore it is useful to calculate J n (q). The result for the 
three different contributions from intra-and inter-band transitions to the susceptibility 
are shown in Fig.([5la-c). The intra-band contributions have their maximum at the Y 
point where the electron part (a) has an additional side maximum at the (3 — (3 nesting 
vector Q = (7r, tt) which is due to excitations between different electron pockets. The 
inter-band contribution (c) on the other hand has maxima at the Q a = (7r, 0), (0,7r), 
i.e. a — (5 nesting vectors. However the value at the T-point (0, 0) is still considerably 
enhanced. For the general FS of Fig. ([He) these contributions have to be summed up 
according to Eq. f |T9l ). The resulting total x*(q) = «/n(q)/ Jl x is shown in Fig. (Eld). Ob- 
viously the absolute maximum is still at the Y point and strongly peaked side maxima 
at Q a = (tTj 0), (0, 7r) are present. This raises a question about the validity of the sim- 
ple parabolic electron-hole model discussed so far since the magnetic instability in the 
3d system of Fe pnictides is of the SDW type with an ordering vector Q a = (7r, 0), (0, tt). 

To understand this issue better we also calculated the momentum dependence of 
Xt(q) — Jn{q)/ Jlx in a more realistic tight binding (TB) type band structure for the 
3d bands valid for the 122 compounds. The Fermi velocities and size of the pockets in 
this model are based on Refs. [30l I3T1 [32] . The band structure and associated density 
of states (DOS) are presented in Fig. ([6]). The detailed parametrization can be found in 
Ref . [47j . The basic features of momentum dependence is similar to the previous model: 
The maxima occur at or close to zone center due to intra-band processes and at zone 
boundary points due to nesting features of electron-electron (a) and electron-hole (c) 
excitations. An essential difference to the parabolic pocket model in Fig.([Hb) is the 
small value of the electron-hole contribution in Fig. ([7b) for small momentum transfer. 
This difference is due to the deep depression in the tight-binding model DOS for uj « 
whereas the DOS for the 2D parabolic band model is simply constant in that region. 
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As a consequence the absolute maximum of Xt(q) = J n (q)/ J 2 X in the tight-binding case 
is located at the Q a = (V, 0), (0, tt) positions ( FiglZld). Therefore in the latter case the 
SDW instability of the itinerant 3d electrons is predicted at the proper wave vector. 

However the simplified model results in Fig.Q are nevertheless useful for the 
understanding of the RKKY interaction because the latter refer to the ordering of 
localized 4f electrons which takes place only at temperatures much lower than the SDW 
transition temperature. In the low temperature region the fully developed SDW gap of 
3d band electrons will strongly suppress the peak in Xt(q) = J n (q)/J 2 x at Q a = (7r, 0) 
(or (0,7r)) because the 'feedback effect' of the gap opening modifies the quasiparticle 
energies connected by these electron-hole nesting vectors and in fact destroys or reduces 
the nesting properties for the quasiparticle bands. Therefore deep in the SDW phase 
the RKKY function will more resemble the one in Fig. (Eld) with the maximum at the 
T point. 

This scenario seems indeed to apply to the 4f- based Fe pnictide compound 
EuFe 2 As 2 . There the AF ordering of itinerant 3d moments takes place below rgjyyy = 
190K at Q a = (V, 0,0), (0,7r, 0) wave vectors [5] It is followed by the AF ordering of 
localized Eu 2+ (S = |) 4f- moments at a much lower temperature T N « 20K and at a 
different wave vector q = (0, 0, 1) which means the Eu 2+ ab-planes with 4f moments 
are ferromagnetically ordered in contrast to the columnar AF order of Fe itinerant 3d 
moments. This agrees with the arguments given above that the RKKY interaction 
within the SDW phase is dominated by the broad peak near the T -point as presented 
in Fig. (Eld). This can be fully understood by the feedback effect on the q-dependence 
of the RKKY interaction in the presence of the SDW ordering (Eq.(l22j)). The SDW 
couples electron (centered at Qi) and hole pockets together which leads to reconstructed 
quasiparticle bands (Sec i3.3l) . Their associated FS has new small pockets around the 
T-point shown in Fig.(Pd). This causes the RKKY contributions of the involved pockets 
and the maximum to move to the zone center. The feedback effect of the SDW not only 
shifts the RKKY maximum to the zone center but also leads to an effective induced 
RKKY spin space anisotropy as argued in Ref. [43] and experimentally supported in 
Refs. [TH HE]. This is also observed in the q-dependence of RKKY as presented in 
Fig. (0a) and Fig. ([Bib) for J x and J z respectively. 

5. Summary 

We have investigated the RKKY interaction mechanism in 3d multiband Fe pnictide 
compounds which is relevant for effective coupling of localized, e.g., 4f moments. We 
used a simple parabolic band model which allows to describe pure electron, hole and 
composite Fermi surface models. In the former cases a closed analytical solution for 
the effective RKKY function was obtained and it is in full agreement with numerical 
calulation. In the latter case with both electron and hole sheets only numerical 
evaluation is possible. The RKKY interaction is determined by the sum of two intra- 
band (e-e and h-h) and two inter-band (e-e and e-h) contributions. The former lead to 
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(a) (b) 

Figure 8. (Color online) The momentum dependence of the Fourier transform of the 
total RKKY interaction in the SDW phase, correspond to a typical tight-binding type 
band structure example for 122 compounds, (see FigEJ): a) J x (q), and b) J z (q). 



slow spatial oscillations determined by the pocket size, the latter to superposed rapid 
oscillations determined by the inter-pocket nesting vectors. Depending on which Fermi 
surface sheets are present only the former or both type of oscillations are present in the 
total RKKY interaction (Fig. [3j). 

We also studied the momentum dependence of the RKKY interaction which would 
determine the magnetic order in 4f-based Fe pnictides in a purely 2D picture. We found 
that in the parabolic band model the maximum of the interaction is at the zone center 
leading to ferromagnetic order of rare earth planes in these compounds in agreement 
with observation. A more closer investigation shows that the RKKY interaction derived 
from a more realistic tight-binding band model would predict a columnar AF order in 
the plane as for the itinerant Fe moments. However, the feedback mechanism of Fe SDW 
gap opening on the RKKY interaction shifts the ordering vector to the zone center as 
obtained in the simplified band model. 

It is worthwile to discuss to what extent this purely 2D model analysis may be 
relevant for the ordering in the 4f based iron pnictides like RFe2As2 and RFeAsO. In 
particular in EuFe 2 As 2 [5] and CeFeAsO the 4f magnetic order has been investigated. 
The former is most directly related to the present model because this 122 structure has 
concerning the local moment properties planar Eu layers. Furthermore the Eu 2+ state 
has a pure spin S = 7/2 without complications due to crystalline electric field effects and 
therefore should closely correspond to the effective spin exchange model. In fact this 
compound exhibits the usual columnar AF SDW order for Fe moments while a much 
lower T N =20 K the Eu- moments in each tetragonal ab plane order ferromagnetically 
as predicted by the present effective exchange model and the feedback effect. The 
stacking along c which is not described in our context is still antiferromagnetic but 
may also become ferromagnetic by substitution of As with P [19J. The 1111- structure 
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corresponds less ideally to our model because the Ce atoms do not reside in planar 
layers. The structure suggests an increased importance of Fe-Ce interlayer coupling 
which would tend to enforce the columnar AF structure also on the Ce layers. In fact 
there are large polarization effects of Ce-moments above the Ce ordering temperature 
due to this coupling [7] . On the other hand the feedback effect still favors ferromagnetic 
effective exchange within the Ce layers. The resulting non-collinear Ce-ordering below 
TV =4 K seems to be a compromise between these two effects; therefore the moments 
of ferromagnetic Ce columns along b are not opposite along a as in the Fe layers but 
only perpendicular. We conclude that at least in these two examples the tendency to 
ferromagnetic effective in-plane exchange found in our model seems to be relevant for 
the real ordering of 4f moments. 
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